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Abstract 

It has recently been proposed that the miscibility of nanoparticles with a polymer matrix 
can be controlled by grafting polymer chains to the nanoparticle surface. As a first step to 
study this situation, we have used molecular dynamics simulations on a single nanoparticle 
of radius R (4a <R< 16a, where a is the diameter of a polymer monomer) grafted with 
chains of length 500 in a polymer melt of chains of length 1000. The grafting density £ 
was varied between 0.04-0.32 chains/a 2 . To facilitate equilibration a Monte Carlo double- 
bridging algorithm is applied - new bonds are formed across a pair of chains, creating two new 
chains each substantially different from the original. For the long brush chains studied here, 
the structure of the brush assumes its large particle limit even for R as small as 8a, which is 
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consistent with recent experimental findings. We study autophobic dewetting of the melt from 
the brush as a function of increasing E. Even these long brush and matrix chains of length 6 
and 12 N e , respectively, (the entanglement length is N e ~ 85) give somewhat ambiguous results 
for the interfacial width, showing that studies of two or more nanoparticles are necessary to 
properly understand these miscibility issues. Entanglement between the brush and melt chains 
were identified using the primitive path analysis. We find that the number of entanglements 
between the brush and melt chains scale simply with the product of the local monomer densities 
of brush and melt chains. 



1. Introduction 

Controlling the interactions, and hence the spatial dispersion, of nanoparticles (NP) is critical 
to the ultimate goal of producing polymer nanocomposites with desired macroscale properties. l ' 2 
Experimental studies have shown that, most often, NPs tend to randomly aggregate into clusters 
or migrate to interfacial regions; consequently, the advantages of their nanoscale dimensions are 
lost in terms of property improvements. Thus, the current barrier to integrating NPs into a range of 
advanced devices is the challenge of controllably dispersing and organizing them within durable 
matrices while retaining their unique properties. One way to alleviate this issue is by end-grafting 
polymer chains to the NPs, a system which is commonly referred to as a polymer brush. 2 For high 
enough coverage, the NPs are sterically stabilized and hence good nanoparticle spatial dispersion 
can result. 5 

Previous work on planar polymer brushes in a chemically-identical matrix has shown that the 
matrix only wets the polymer brush if the melt chains are shorter than the brush. 611 Longer melt 
chains are known to spontaneously dewet the brush ("autophobic dewetting"), because the entropic 
loss associated with the matrix chains penetrating the brush cannot overcome the translational, or 
mixing, entropy. 5 ' 12 ' 13 Thus, for a given graft density, increasing the chain length of the polymer 
matrix effectively changes the solvent quality experienced by the polymer brush from "good" to 
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"poor". Similar ideas il ' 14 ~ lb can be used to understand the miscibility of polymer grafted NPs in a 
polymer melt, except that, in addition to the graft density and chain lengths of the brush and melt, 
one must also account for the curvature of the NP. For small NPs (high curvature), the polymer 
brush chains can explore more space compared to a planar surface, resulting in less entropic loss to 
penetrate the brush, reducing the tendency for autophobic dewetting. Most previous computational 
studies of polymer brushes grafted to a nanoparticle have either not included the polymer matrix 2 ' 17 
or have studied only short chains in which both brush and melt are well below the entanglement 
chain length, A^. 18 There have also been studies on bare spherical nanoparticles in long chain 
polymer matrices. 19 

Perhaps the most commonly accessible measure of brush structure is its mean height, a quantity 
that has been the focus of several experimental studies. Chevigny et al. 20 recently used neutron 
scattering to measure the conformation of silica NPs with an average radius of 13 nm grafted with 
a polystyrene polymer brush dispersed in a polystyrene matrix. They found that if the melt chains 
become long enough, the brushes compressed by a factor of two in thickness compared to their 
stretched conformation in solution. Other experiments, 21 " 23 which also measure the mean brush 
height, have been able to critically comment on existing theories of curved polymer brushes. For 
example, theories that assume that all the free chain ends are a uniform distance from the grafting 
surface do not provide a good representation of the experimental results. 24 ' 25 Theories which relax 
this assumption 14-16 provide a much better description of the experimental data. These theories 
also predict that there is a exclusion zone for the free ends - that is, the ends cannot approach 
the particle surface, thus creating an "exclusion zone". This result is in sharp contrast to planar 
brushes, where no exclusion zone is predicted from theory or found in experiments. 

While previous work has focused on a single moment of the spatial distribution of monomers, 
the mean brush height (the first moment), examining the full distribution of monomers is a more 
sensitive and critical test of the theories. This type of comparison is best performed via simulation 
studies, which we undertake here. To provide context, we consider the theoretical predictions of 
radial distance dependence of the density of grafted chain monomers. 
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In the poor-solvent limit, Daoud and Cotton predicted that the brush chains collapse into a 
globule whose chain density profile is simply a constant value of where v and ft) are the 
second and third virial coefficients respectively. 24 This result is independent of the curvature of 
the particles and reflects the dominance of the solvent conditions alone. For © and good solvents, 
in the limit of large curvature (small particle radius, R), the chain density profile is predicted to 
scale as 15 



where r is the radial distance from the center of the particle, E is the graft density and the exponent, 
m, is 1 or 4/3 for a or good solvent, respectively. Note that is nonlinear in E because the brush 
height H depends on E. Implicit in this treatment is the unrealistic assumption that all chain ends 
are extended to an identical degree, and fixed to a "phantom surface"; an assumption that is not 
representative of physical systems. In contrast, in the limit of small curvature (large R) one must 
account for a distribution of chain ends, e.g., as proposed by Milner-Witten-Cates. 26 Following 
this form, the segment density profile is parabolic, given by 15 



where a is the statistical segment length, z = r — R is the radial distance from the particle surface, 
H is the effective brush height for a flat surface, and H is the brush thickness. Note that this form 
is reminiscent of planar polymer brushes and suggests that, asymptotically, we must recover the 
planar brush result for arbitrarily large particles. 

Finally, for intermediate particle sizes, a combination of large and small curvature behaviors 
is conjectured. Thus, the segment density profile has large curvature behavior near the surface 
[Eq. (1)], followed by small curvature behavior away from the surface. The crossover between the 
two behavior occurs at zo, which is obtained by matching the two solutions. The combined profile 
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where / [= L(4nR 2 )] is the number of grafted chains, v is the second virial coefficient, z is the 
location of the transition, N2 is the number of segments beyond z , and H2 is the brush thickness 
beyond z - 

In addition to the form of the monomer distribution function, these theories also predict the 
density profiles of the (free) chain ends. The Daoud-Cotton model assumes that all the chain 
ends are a uniform distance away, implying a delta function distribution. The Wijmans-Zhulina 
model, the spherical brush analog of the Milner-Witten-Cates model, has a well defined exclusion 
zone (where the probability of finding a chain end becomes negative and hence unphysical). The 
expression for the density of free ends is given by eq. 19 in ref. 15 

For the long chains modeled in this study, we can readily determine if entanglements occur 
between the brush chains and the matrix chains. Such information is useful in determining the 
role of these entanglements in the mechanical reinforcement gained by a grafted particle immersed 
in a polymer melt. A particularly useful concept in this context is that of the primitive path, 27 
which gives a measure of how the presence of other (long) chains forces the central chain to be 
transiently confined in a "tube". Using the primitive path analysis (PPA) developed by Everaers and 
co-workers, 28 it is now possible to get an approximate measure of an individual chain's primitive 
path, entanglements, and hence interactions with all other chains in the melt via relaxation methods 
of individual melt configurations. 

In this work, we use molecular dynamics (MD) simulations to model spherical nanoparticles 
for three values of radii, R, and four coverages, E, coated with end-tethered chains of length N t , 
immersed in a melt of chains of length N m , where both N t and N m are much greater than the entan- 
glement length N e . In the next section, we present the polymer model and simulation methodology. 
We then present results for the monomer density profiles for both the brush and melt chains. The 
conformations of the brush and melt chain are examined at various scales and compared to theory. 
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Finally, using PPA analysis, 28 ' 29 we examine the properties of the system entanglements within 
the polymer brush and between the polymer brushes and the melt. 

2. Model and Simulation Methodology 

2.1 Interaction potentials and system parameters. End grafted polymer coated NPs in a poly- 
meric matrix are simulated using the coarse grained bead spring model of Kremer and Grest. 30 All 
brush and melt monomers have a mass m and interact via the truncated and shifted Lennard- Jones 
potential 



r > r c 

where o it the monomer diameter and e is the characteristic pair interaction energy. Covalently 
bonded monomers also interact via the finitely extensible nonlinear elastic (FENE) potential 

kR 2 

U FE NE = —fln{l-(r/R ) 2 ), (5) 

where k = 30ksT and R Q = 1.5c 7 The NP is modeled as a spherical particle of radius R consisting 
of a uniform distribution of monomers of diameter o and density p NP . The interaction between 
the NP and the brush and melt monomers is given by integrating Eq. (4) over this sphere: 31 
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where Anp, p = 24n£ nv p N p(7 3 . Here £ np = £ is the interaction between a polymer monomer and a 
monomer in the nanoparticle. p NP = 1 .0a~ 3 gives A NP n = 75.4c. The NP-monomer interaction is 
truncated at r c = (R + a). 

For all systems, n t chains were randomly tethered at a coverage E to the surface of the NP, with 
the only restriction that no two grafting points lie within 2a of each other. Three values of the NP 
radius were studied, R = 4, 8, and 16c, for coverages E = 0.04, 0.08, 0. 16 and 0.32 chains/cr 2 . The 
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number of tethered chains varied from 8 to 1024 depending on the coverage and NP radius. The 
number of melt chains n m = 2000 for R = 4a, 3000 for R = 8 a and 5000 for R = 16(7. Each melt 
chain has a degree of polymerization N m = 1000, while each tethered chain is of length N t = 501 . 
For flexible bead spring chains, the entanglement length N e for bulk melts is N e = 85 ±7, 32 so the 
polymer melt is well entangled. 

For NP-brush-melt systems, the initial states were created as described in Auhl et al. 33 with 
the additional restriction that no monomers overlap with the NP which is held fixed at the center of 
a periodic, cubic LxLxL simulation cell. L was initially chosen so that overall monomer density 
p m = 0.75<7~ 3 . Overlapping monomers in the initial states were pushed off each other using 
a soft potential until the monomers were far enough apart to switch on the full Lennard- Jones 
interaction. After the full Lennard- Jones potential was turned on the system was equilibrated at a 
pressure P = 5.0e/a 3 which for this system is the pressure of a neat polymer melt of long chains 
at a monomer density p m = 0.85a -3 . 30 The equilibrated size of the simulation cell is L ~ 132a 
for R = 4a, ~ 154 for R = So and ~ 182a for R = 16a, with slight variations depending on the 
coverage E. The number of melt chains was sufficiently large to ensure that the tethered chains do 
not interact with their periodic images. 

The systems in this study were thoroughly equilibrated using the double-bridging method de- 
scribed in Refs. 29 ' 33 " 35 Because diffusive equilibration of entangled brush-melt systems is expo- 
nentially slow, Monte Carlo moves that alter the topological connectivity of chain subsections were 
periodically performed, allowing the chain configurations to relax more rapidly. 33 ' 34 

We also simulated NP-brushes in an implicit solvent, as has been done in previous simulations 
of a dry brush on a flat surface, 36 ' 37 a cylindrical surface 38 and a sphere. 17 All simulations of the 
NP-brush-melt systems and NP good-solvent brushes are run at T = e/kg and r c = 2 1 / 6 a. We also 
simulated NP-brushes in an implicit poor solvent with T = 1 .5£/fcg and r c = 2.5a. For this model, 
the © temperature is T© ~ 3.18£/fcg. 39 

All simulations were carried out using the LAMMPS parallel molecular dynamics package. 40 
Newton's equations of motion were integrated with a velocity- Verlet algorithm with a time step 
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8t = 0.012t for the melt and good solvent simulations and 0.005T for the poor solvent, where 
T = a(m/g) 1 / 2 . All monomers were coupled to a Langevin thermostat with a damping constant 
r = 0.5t _1 . Figure 1 shows representative snapshots from the simulations, where it is clear that 
there is intimate contact between the brush chains and the surrounding melt. Increasing the graft 
density serves to sterically shield the NP surface, suggesting that brush physics should ultimately 
decide the thermodynamic state (miscible vs. clustered particles) of these composites. 




Figure 1: Snapshot of a NP of radius R = Ao for (a) E = 0.04 and (b) E = 0.32 chains/a 2 . Most 
of the melt chains (blue) have been removed to visualize the grafted chains (red). 

2.2 Primitive Path Analysis For NP-brush-melt systems, we performed a primitive path analysis 
(PPA) to identify the entanglements between the brush chains and between the brush and melt 
chains. 28 ' 29 All chain ends are fixed in space, and several changes are made to the interaction po- 
tential. Intrachain excluded- volume interactions are set to zero, while interchain excluded- volume 
interactions are retained. The covalent bonds are strengthened by setting k = lOOe/d 2 . FENE 
bonds are monitored to ensure they do not exceed 1.2a; this prevents chains from crossing. 41 Self- 
entanglements are not necessarily preserved, but their contribution is typically negligible. 41 ^ 3 The 
system is coupled to a heat bath at T = 0.001 <J /kg, and the equations of motion are integrated until 
the chain lengths are minimized. 

Since the standard PPA is not sufficient to identify individual entanglements, following the 
convergence of the standard PPA, we therefore perform "thin-chain PPA" using the same method 
described in Ref . 29 At the end of the standard PPA, n dec — 1 new beads are inserted between 

8 




Figure 2: Schematic of entanglement identification. For visualization purposes, chain A is red, 
chain B is blue, and chain C is green. 

adjacent beads on the original chains. Here, n^ ec = 8. The bead diameter (i.e., range of the LJ 
potential), range of Lennard Jones interaction between the nanoparticle and system monomers, 
and Rq are reduced to l/n^ec of their initial values. Covalent bonds are further strengthened by 
setting k = \00n c [ ec /o 2 , and their lengths are capped at 1.2cr/n^ ec to avoid chain crossings. The 
equations of motion are integrated until the chains again minimize their length. Results presented 
here are averaged over 10 to 20 statistically independent systems. 

Figure 2 schematically depicts our PPA analysis procedure. One entanglement between chains 
A and B is identified as a block of consecutive monomers on chain A having interchain contacts 
with chain B. As described in Ref. 29 identifying these blocks is somewhat subtle; we identify 
interchain "contacts" as monomers whose smallest interparticle neighbor distance is less than 
fcon = (5/4)(l/n^ ec )2 1 / 6 a. For computational simplicity, we follow convention and assume all 
entanglements are binary, i.e., involve only two chains. Ternary entanglements (a vertex including 
three chains A, B, C) are counted as three entanglements: A — B, A — C, and B — C [Figure 2]; 
this may overestimate the absolute entanglement densities, but the relative values of brush-brush 
versus brush-melt entanglement densities should be accurately captured. 
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3. Chain Structure 

The radial density, (j>(r), of the tethered chains (Figure 3 a-c) and the corresponding melts 
(Figure 3 d-f) are shown for three different nanoparticle radii R = 4, 8, and 16a, respectively, 
and four different grafting densities, E = 0.04, 0.08, 0.16, 0.32 chains/a 2 . It is apparent that the 
radial monomer density profile of the brushes extends farther into the melt with higher grafting 
densities. Further, the shape of the monomer density profile changes quite dramatically as the NP 
radius increases from R = 4 to R = 16c The smallest NPs show a concentration profile that is 
more concave, while for the largest NP the profiles are more parabolic in its shape, especially for 
intermediate values of r. 
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Figure 3: Monomer density profiles for NP radius R = 4a (parts a and d), 8c (b and e) and 16c 
(c and f). The density of grafted chains is shown in the top row while the density of melt chains is 
shown in the bottom row. 



The melt density profile (Figure 3d-f) typically mirrors the density properties of the grafted 
brushes. The sum of the two densities suggests that the system is essentially incompressible, 



10 



except in the immediate vicinity of the grafting surface. This result is not surprising given the high 
overall density of the system. 
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Figure 4: Comparison of the normalized segment density profiles to large curvature behavior for 
a) R = 4a, b) R — 8 a, and c) R = 16a. Slopes of 1 and 4/3 (cf. Eq. (1)) are drawn as a guide to 
the eye. 

To understand the dependence of the monomer density profile on the E and R, we perform an 
analysis inspired by Daoud-Cotton, i.e., following Eq. (1). Figure 4 shows that this scaling only 
approximately describes the density profile at small r or large R/r. For all particle diameters and 
small E it appears that a scaling exponent of m = 4/3 approximately describes the data near the 
particle surface. As E increases, the solvent quality appears to get progressively poorer, evidenced 
by the decay in the scaling exponent. This suggests a transition from "good" to "poor" solvent con- 
ditions with increasing grafting density. Note, however, that these ideas are only qualitative since 
there is no region where the scaling relationship is followed exactly; rather, the data are always 
curved. For large r (small R/r), however, the data do appear to be better fit by an approximately 
exponential behavior reminiscent of that found for flat surfaces - this behavior presumably follows 
from the finite chain lengths, and hence finite extensibility, of the graft chains simulated. We can 
thus conclude that the model of Daoud and Cotton only qualitatively describes the simulation re- 
sults for the brush density profile. However, it evidently shows that solvent quality becomes worse 
as one considers more densely grafted particles, a situation in which one may expect a smaller 
propensity for the melt chains to wet the brush. 
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Figure 5: Comparison of the normalized segment density profiles to small curvature behavior ( 
Eq. (2)) for a) R = 4<T, B) R = So, and C) R = I60. 

To obtain a better understanding of the brush density profile, we also examined the applicability 
of the large particle (large R) ideas embodied in Eq. (2). Figure 5 shows that Eq. (2) surprisingly 
describes the results obtained for all R except for R = 4a. While this treatment describes the 
behavior for all small r values for R > So, apparently for all E values we consider, poor agreement 
is found for larger r values. As described above, this discrepancy probably results from the finite 
extensibility of the chains, a factor not built into these theories. Consistent with these findings 
Figure 6 shows that the transitional form Eq. (3), which is supposed to apply for intermediate R 
values, does not provide improved fits. An important conclusion we draw from this analysis is that 
the large particle limit is attained relatively quickly, and particles as small as 10 times the monomer 
diameter of the chains are already in this limit. 

A few experimental points are considered in the context of these simulation results. One ques- 
tion: is the apparent good-solvent behavior found for the largest R reasonable, given that the brush 
should "collapse" due to the autophobic dewetting that is expected for long enough matrices and 
high enough grafting densities? Note here that all of the monomer density data for R > So appear 
to be fit by a parabolic profile which is indicative of good solvent behavior. Thus, at first glance 
there appears to be no signature of the dewetting phenomena in these monomer density profiles. 
We then conjectured that, perhaps, a moment of the monomer density, namely the mean height 
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Figure 6: Comparison of the normalized segment density profiles to intermediate curvature behav- 
ior (Eq. (3)) for a) R = 4a, b) R = 8a, and c) R = 16a. 
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might give better insights into this dewetting phenomenon (Figure 7). In addition to this "di- 
rect" measure of height, another measure can be derived by fitting the monomer density profile 
to well known functional forms. In particular, we follow an approach that we previously used 
to analyze the effective % parameter between the brush chains and the matrix. 13 We first define 
a normalized density, p = ^ gra ft(z) / (^ gra ft(z) + ^matrix(z))- In ref. 13 we showed that concen- 
tration profiles derived from self-consistent mean-field theory could be fit to a standard form: 



1 — tank 
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where z is the position of the interface and W\/2 is the interfacial width. 
Paralleling the approach for phase segregated diblocks, then, W\n = (6#)~ 1//2 . This form fits all 
of our data very well, and in Figure 7 we report both the z and W\/2 values derived from the fits. 
Clearly z tracks the behavior seen for the mean brush height < z >. The results (Figure 7a) show 
that for the largest radius particle, < z >~ E 1 / 2 , is consistent with the expectation of a stretched 
brush in a good solvent. In contrast, the R = 4a data are very noisy, but are essentially indepen- 
dent of E. More interestingly, the interface width W\/2, in all cases, converges to a value of around 
10- 12a at the largest surface coverages, E, examined. This number is comparable to the radius 
of gyration of the melt chains, R g ~ a/ olKN m /6 ~ 17a (where the Kuhn length Ik ~ 1.8a), sug- 
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gesting that we are limited by the size of the chains in question. An unambiguous measurement of 
interfacial width, without interference from this length scale, requires that the two length scales be 
very different. Even though the chains employed here are relatively long, they are still too short 
to allow for this separation of length scales and hence we cannot determine if the brush and melt 
chains are phase- separating due to autophobic dewetting. 

In this context, it is also informative to compare the melt systems with simulation in an implicit 
good and poor solvent. For a typical system (see Figure 8), close to the surface, we see that the 
curved brushes interacting with a melt have a density profile that resembles more that seen for an 
implicit, poor solvent than for a good solvent. However, there is a much longer tail for the brushes 
immersed in a melt. Similar behavior was also previously observed for brushes grafted to a plane; 
in both cases the 'long tail' is probably attributable to finite chain length effects. 

A particularly interesting prediction of mean-field theories of curved brushes is the existence 
of a zone of exclusion for the free chain ends near the surface. To study this interesting possibility 
we examine the free chain end distributions, p e (z) of the polymer brush in Figure 9. It is apparent 
that there is a maximum in this quantity, but that there is no complete exclusion of chain ends 
from near the particle surfaces. These results can be compared to the predictions of the Wijmans- 
Zhulina theory in the limit of large particles. The theory does however captures the maximum in 
the end density profile at some distance away from the particle surface and its relative position with 
increasing graft density. 

Finally we identified the entanglements between the polymer brush and the melt. Figure 10 
shows results for the brush-brush, brush-matrix and matrix-matrix entanglement densities, Pbb(z), 
Pbm(z) and Pmm(z) respectively scaled by the melt density p m for E = 0.04 — 0.32 chains/a 2 and 
R = 4o. The densities Pes(r) and pBM(r) increase with increasing coverage E. The maximum 
entanglement density, which occurs at the same r where the total monomer density also has a peak, 
scales as E L4 . As previously observed for planar brushes, 29 the fraction of brush-melt entangle- 
ments scales to a good approximation as the product of the monomer densities, §brush{ r )§melt{ r )i 
as expected for 'binary' entanglements. Figure 10b shows that the fraction of brush-melt entan- 
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Figure 7: (a) Mean brush height < z>, (b) interface position z , and (c) interface width W\/2 as a 
function of coverage E for NPs of radius R = 4 (circles), 8 (squares), and 16c (triangles). 
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Figure 8: Brush monomer density for NP of radius R = 8<7 for (a) E = 0.04 chains/a 2 and (b) 0.32 
chains/a 2 in a polymer melt and in an implicit good and poor solvent. 



16 



0.0025 



17 0.002 

| 0.0015 

<D 
Q 
"D 

LU 
C 

O 0.0005 



1 r 

a) R=4o 



1 1 1 r 

b) R=8o 



0.001 



ilMi 



• 


0.04 


♦ 


0.08 


■ 


0.16 


▲ 


0.32 






10 20 30 40 50 60 



10 20 30 40 50 60 10 20 30 40 50 60 

z z 



Figure 9: End density p e (z) of the tethered chains in the presence of melt chains for NP radius (a) 
R = 4a, (b) 8a, and (c) 16c for four values of the coverage E. The lines are predictions of the 
Wijmans-Zhulina theory. 15 



glements peaks at a distance 3 — 10c away from the surface of the nanoparticles. The z position 
of the maximum increases with increasing E. Similarly, the fraction of melt-melt entanglements 
(see Figure 10c) also scales roughly quadratically with the melt monomer density, and increases 
monotonically with increasing z. Given the simple behavior observed, we did not consider the 
larger radii where similar results are expected. Results for for planar brushes are reported in Ref. 29 
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Figure 10: Scaled densities for (a) brush-brush entanglements, (b) brush-melt entanglements, and 
(c) melt-melt entanglements for a NP of radius R = 4a. All distances are measured from the NP 
surface. Results are averaged over 10-20 independent initial configurations. 
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4.0 Conclusions 

We have performed detailed MD simulations of the behavior of a polymer brush grafted to a 
single NP in contact with a long chain melt. Our results have a number of salient points, (i) We 
show that the brush monomer density for NPs of radius R = Ac is close to the behavior anticipated 
for a star polymer, while, for radii R = 8 and 16c the brush monomer density is already close to 
that of the large particle limit. This result is consistent with experimental results, (ii) The con- 
formational transition between a stretched brush and a collapsed brush, which has been attributed 
to the autophobic dewetting of the brush by the melt, is not readily apparent in the brush density 
profiles. We conjecture that this difficulty is due to the finite size of the chains in question, even 
though the number of beads in the simulated brush chains is long enough to be well entangled - 
we, thus, cannot cleanly resolve the interfacial width attributed to dewetting effects for the native 
sizes (R g ) of the chains in question. Simulations of two or more grafted nanoparticles, which will 
help to unequivocally resolve this issue, are currently underway, (iii) The brush end-monomer 
density profiles are qualitatively consistent with the large nanoparticle predictions of Zhulina and 
Wijmans, though we do not see any clear evidence for a zone of depletion. This may also be a 
manifestation of the finite chain lengths considered in our work, (iv) The spatially dependent en- 
tanglement densities semiquantitatively scale as a product of the respective densities of the two 
chains that form the brush-brush, brush-melt and melt-melt entanglements respectively. 
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